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Abstract 

We study the time evolution of freely rolling moduli in the context of M-theory on a Gi manifold. 
This free evolution approximates the correct dynamics of the system at sufficiently large values 
of the moduli when effects from non-perturbative potentials and flux are negligible. Moduli fall 
into two classes, namely bulk moduli and blow-up moduli. We obtain a number of non-trivial 
solutions for the time-evolution of these moduli. As a generic feature, we find the blow-up moduli 
always expand asymptotically at early and late time. 
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1 Introduction 



Seven-dimensional manifolds with holonomy G2 provide the general framework for relating M-theory 
to four-dimensional models with N = 1 supersymmetry. Moreover, it has been shown that for 
certain singular limits of such manifolds four-dimensional models with phenomenologically interesting 
properties can be obtained [1]— [16]. 

Given these promising features it seems worthwhile to have a closer look at the cosmology of such 
G2 compactifications of M-theory. More specifically, in this paper, we will be interested in the time- 
evolution of G2 moduli. This problem has not previously been analysed since the four-dimensional 
effective theories, specifically the kinetic terms, for these fields were unknown. However, in a recent 
paper [17], the four-dimensional moduli Kahler potential for M-theory on a G2 manifold has been 
computed explicitly and our analysis will be based on these results. 

The G2 manifold considered in Ref. [17] is based on the compact examples due to Joyce [18]-[20], 
constructed by blowing-up a seven-dimensional G2 orbifold. The moduli for such a G2 manifold 
split into two groups, namely the bulk moduli (related to the radii of the underlying seven-torus) 
and the blow-up moduli (related to the volume of the blow-ups). The Kahler potential of Ref. [17] 
was obtained in the usual large-radius supergravity approximation and for blow-up volumes small 
compared to the total volume of the space. In the limit of a vanishing blow-up modulus the space 
develops a singularity (which corresponds to a singularity of the underlying orbifold) leading to 
new fields in the low-energy theory. It is singular limits such as this which are essential for a 
phenomenologically interesting low-energy theory. 

In this paper, we will focus on the simple case of free moduli evolution based on the Kahler 
potential of Ref. [17]. That is we will neglect the effect of possible non-perturbative potentials and 
flux. Hence our results constitute a good approximation at large moduli values where effects of 
a non-trivial moduli (super)-potential can be neglected. Particular consideration is given to the 
evolution of the blow-up moduli; particularly as to whether they may contract, thereby approaching 
one of the phenomenologically interesting singular limits. 

Let us summarise our main results. Starting from the given Kahler potential and consistently 
truncating off the axions in the chiral moduli multiplets, we show that, after a suitable reparametri- 
sation of fields, the sigma-model metric becomes independent of the blow-up moduli. As a result, 
these fields can be integrated out and the problem reduces to one of studying the evolution of the 
bulk moduli in an effective potential. Based on this approach, we show that a number of explicit 
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analytic solutions can be obtained using methods from Toda theory. For all these solutions, we find 
the relative volume of the blow-up moduli increases asymptotically at early and late times. As a 
consequence, the small blow-up approximation underlying the Kahler potential breaks down after 
a finite proper time. This feature can also be intuitively understood from the properties of the 
effective bulk-moduli potential and can, therefore, be viewed as generic. Hence, once moduli are 
large it appears to be difficult to evolve towards a singular state of the internal space. This suggests 
that a successful cosmological evolution leading to phenomenologically interesting singular G2 spaces 
should stabilise moduli at small values from the outset. 



2 Kahler potential for G2 moduli 

In Ref. [17] the moduli Kahler potential has been calculated for a specific manifold with G2 holonomy 
constructed from a G2 orbifold by blowing up the orbifold singularities. The chiral moduli multiplets 
for this manifold split into two categories, namely the bulk moduli T A , where A, B,C, ■ ■ ■ = 1, . . . , 7, 
which encode the radii of the seven-torus T underlying the orbifold, and the blow-up moduli U l , 
which measure the volume of the blow-ups. The general structure of this Kahler potential is 

7 

K = - ^2ln(T A + f A ) + ^2 Fi{T A ,f A ) (U l + U i ) 2 + c , (2.1) 

A=l i 

where the functions Fj are given by 

and c is a constant. Here, Ai and Bi specify the two particular bulk moduli by which the blow-up 
modulus i is divided in the above Kahler potential. The values for these indices reflect the underlying 
structure of the orbifold and they can be conveniently encoded in constant seven-dimensional vectors 
Pi by writing 

^ = 2e- pf3 (2.3) 
where the new fields (3 = ([3a) are defined by 1 

Re(T A ) = e pA . (2.4) 

We expect this form of the Kahler potential to apply to a wide class of G2 manifolds constructed by 
blowing up orbifolds, with the specifics of each example encoded in the vectors p^. 



1 Here and in the following bold-face symbols denote seven-dimensional vectors. 
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P(t,o) 


a = 1 


a = 2 


a = 3 


t = a 


(1,0,0,0,0,1,0) 


(0,1,0,0,1,0,0) 


(0,0,1,1,0,0,0) 


T = (3 


(1,0,0,0,0,0,1) 


(0,0,1,0,1,0,0) 


(0,1,0,1,0,0,0) 


T = 7 


(1,0,0,1,0,0,0) 


(0,0,1,0,0,1,0) 


(0,1,0,0,0,0,1) 



Table 1: Values of vectors P( T , a ) which define the moduli Kahler potential for the G2 manifold based 
on T 7 /Zl 



In the following we will, for concreteness, focus on the particular example of Ref. [17], which is 
based on the orbifold T 7 The blow-ups for this particular manifold can be labeled by a triple 
(i) = (r, n, a) of indices where r = a, j3, 7 indicates the type, that is, under which of the three Z2 
orbifolding symmetries the associated fixed point remains invariant, n = 1,2,3,4 labels the fixed 
points of equal type and a = 1,2,3 is an index associated to each fixed point that describes the 
orientation of the blow-up relative to the bulk. Hence, there are 12 fixed points, labeled by (r, n), 
each with three associated blow-up moduli and, consequently, 36 blow-up moduli in total. For 
simplicity of notation, we will use the single index i to label the blow-ups whenever possible and 
only split into the triple (r, n, a) when required. It turns out that the vectors p« for this case only 
depend on the type r and the orientation a. The resulting nine vectors are given in Table 1. 

It will be useful for the subsequent discussion to have an interpretation for the moduli and the 
various parts of the Kahler potential (2.1) in terms of the geometry of the underlying G2 manifold. 
To this end, one notes that the moduli Kahler potential for M-theory on a G2 manifold is generally 
related to the volume V of this manifold by [10] 

tf = -31n(^). (2.5) 

The first, T^-dependent, part of the Kahler potential (2.1) then corresponds to the volume Vb of 
the orbifold, while the second, [^-dependent, part measures the reduction of this volume due to the 
blow-ups. For the orbifold volume Vq one easily finds from Eqs. (2.5) and (2.4) that 



V = -exp - Y8 A ] , (2.6) 



where we have used the value c = 61n(87r) + ln(2) found in Ref. [17]. From a similar calculation one 
finds that the fraction by which the blow-up (r, n) reduces the orbifold volume Vq is given by 

Hr,n) = e-*<-)*4 >B>a)> (2.7) 

a=l 



where we have introduced the real parts 



Re(f/ i ) = Ui (2.8) 

of the blow-up moduli. 

The Kahler potential (2.1) relies on two approximations so we should discuss its range of validity 
in moduli space. First, we require all moduli to be larger than one, that is 

T A > 1 , 17* » 1 , (2.9) 

for the supergravity approximation underlying the calculation of K to be valid. Second, the Kahler 
potential has been calculated to leading (quadratic) order in U l /T A and terms of order four in these 
ratios have been neglected. Consequently, application of the Kahler potential (2.1) should be confined 
to the region of moduli space where the ratios U l /T A are smaller than one, or more precisely, where 

Hr,n) « 1 • (2-10) 

From the above interpretation of e( T , n ) this means that the volume taken away by the blow-ups should 
be small compared to the volume of the orbifold. In subsequent calculations, we will consistently 
apply this approximation in that we neglect higher order terms in €( T ,n)- 



3 A simple form of the Lagrangian 

The Kahler potential described in the previous section depends on 43 chiral superfields. We, therefore, 
have 43 real moduli fields from the real parts of theses superfields, which are associated to the 
geometry of the G2 manifold, plus 43 axions from the imaginary parts. It is clear that the axions 
can be consistently set to constants and for simplicity this is what we will do in the following. We 
are, hence, left with the real scalar fields Pa and Uj. 

However, their kinetic terms, as computed from Eq. (2.1), are still fairly complicated. A sub- 
stantial simplification can be achieved by introducing a new set of fields 4>a and Zi defined by the 
Ansatz 

Pa = 0a + 5>Ae qr< ^ 2 t 3 - 1 ) 

i 

m = e 3 '^ . (3.2) 

Here, Cj = ((Ha), q« = (lAi) and Sj = (sAi) are constants to be determined shortly. The kinetic terms 
for Pa and u,- t can now be rewritten using this field reparametrisation. For the particular choice of 
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constants 



Cj = -4/XjPi (3.3) 

= i(^ + l)p 4 (3.4) 

q» = MVi , (3.5) 

where \ii € {—1, +1} are arbitrary signs, we find that the Lagrangian simplifies to 



(3.6) 



where 

P, = . (3.7) 

The obvious advantage of this form is that the rescaled blow-up moduli Zi (unlike their counterparts 
Ui) only appear through their derivatives. This implies that their equations of motion can be imme- 
diately integrated once, leading to a set of first integrals. We will use this property in the following 
section when we discuss the cosmological evolution based on the Lagrangian (3.6). Since we have 
neglected higher order terms in £( Tjn ) in its derivation, it is equivalent to the one directly obtained 
from the Kahler potential (2.1) only provided fields are in the region of moduli space defined by 
Eq. (2.10). Also note that this equivalence holds for an arbitrary choice of the signs /ij. 

4 Cosmological evolution equations 

To study the cosmological evolution based on the Lagrangian (3.6) we start with a metric 

ds 2 = -e 2 ^dr 2 + e 2a ^d^ 2 (4.1) 

of Friedmann-Robertson- Walker type, with the spatial sections taken to be flat for simplicity, and 
time-dependent moduli fields 

4>A = 0a(t) , Zi = Zi ( T ) . (4.2) 

Here, a(r) is the scale factor of the universe and v(t) is the lapse function, which we will determine 
later. From Eq. (3.6) the equations of motion for the blow-up moduli z% are given by 

where the dot denotes the derivative with respect to r. They can be easily integrated to 

Zl = Cie F ^- 3a+u , (4.4) 
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with arbitrary integration constants Q. Let us now choose the particularly convenient gauge v = 3a. 
In this gauge, the sum of the (00) component and the spatial components of the Einstein equations 
reads a = and, we have 

a = v t + Qo (4-5) 
where v and ao are arbitrary integration constants. Proper time t is obtained by integrating 

dt = e 3a dr (4.6) 

which leads to 
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+ a = iln(v^B|t-t |) • (4.7) 
Here to is another integration constant and 

E = 3v 2 . (4.8) 

The solution (4.5) for the scale factor, written in proper time, then takes the form 

a = - hi (y/SE\t - to\j (4.9) 

and shows the expected power-law behaviour, with power 1/3, characteristic of evolution driven by 
kinetic energy. As usual, we have two branches, the negative-time branch, t — to < 0, which ends in 
a future curvature singularity, and the positive-time branch, t — to > 0, which starts out in a past 
curvature singularity. 

We still have to consider the equations of motion for the bulk moduli (f>A- Replacing the blow-up 
moduli with Eq. (4.4) and working in the gauge v = 3a they take the form 

The "effective" potential 

V = 4^(fe F ^ (4.11) 

i 

in these equations of course originates from integrating out the blow-up moduli. The above equations 
of motion can be obtained from the Lagrangian 

£ = ^0-0-y. (4.12) 

They have to be supplemented by the Hamiltonian constraint 

n = -4>-4> + v = e , (4.13) 



which is simply the Friedmann equation rewritten in our language. 



For the subsequent discussion it is useful to express the quantities measuring the geometry of the 
Gi manifold in terms of the redefined fields. The volume Vo of the underlying orbifold, defined in 
Eq. (2.6), can then be written as 




(4.14) 



where we have neglected higher order terms. Likewise, the fraction e( T ri ) of the total volume taken 
up by the blow-up (r, n) now takes the form 

£ (r,n) = |E e_Pi "*4,»,a)- ( 415 ) 
a 

Let us summarise what we have achieved so far. We have integrated out the blow-up moduli 
and decoupled the scale factor by choosing a particular gauge. This has reduced the problem of 
analysing the evolution of 43 freely rolling G2 moduli to one that only involves the seven bulk 
moduli subject to the potential (4.11). The remaining problem is, therefore, to solve the seven 
equations of motion (4.10) for the bulk moduli and the constraint (4.13). Each solution for the bulk 
moduli can then be inserted into Eq. (4.4), which determines the corresponding evolution of the 
blow-up moduli. 

There exists a well-developed solution theory [21], [22]-[26] for the Lagrangian (4.12) and the 
associated equations of motion (4.10) for the bulk moduli if the system is of Toda type. To discuss 
this more specifically we introduce the matrix 

A ij = P t -P j , (4.16) 

which consists of one or more irreducible blocks along the diagonal. The system is called Toda if 
each of these blocks is proportional to the Cartan matrix of a simple Lie-group. One can now inspect 
the vectors P, = —/^iPi from Table 1, where we recall that the fii are arbitrary signs. It is clear that 
the complete nine-dimensional matrix obtained from (4.16) is not proportional to a Cartan matrix 
for any choices of the signs ^ and, hence, the system is not of Toda type in this case. However, 
note that the potential (4.11) does not necessarily contain all possible terms. In fact, each term is 
multiplied by an arbitrary integration constant Q that can be set to zero. Such a vanishing Q implies, 
from Eq. (4.4), that the corresponding blow-up modulus Zi is constant. Hence, we learn that we can 
consistently freeze an arbitrary subset of blow-up moduli Z{ and thereby select an arbitrary subset 
of the vectors P« that appear in the potential (4.11). Then, for suitable subsets and choices of signs 
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Hi, the associated matrix A can well be proportional to a Cartan matrix, as we will see for several 
explicit examples studied further below. The system is then Toda and can be integrated analytically. 
The resulting solutions are, of course, special in that they rely on a number of blow-up moduli being 
frozen. However, in the present context, this seems to be the best one can do by analytic methods. 

Before we go on to analyse explicit examples, it is worth making a general observation about 
the structure of the solutions. It is well-known [24, 25], for Lagrangians of the type (3.6), that the 
fields Zi approach constant values asymptotically at early and late times (both in the negative and 
the positive time branch). This behaviour will also be confirmed in the explicit examples. Further 
note that the potential (4.11) consists of a sum of positive terms. Therefore, one expects the modes 
Pj ■ <j> to roll up the exponential slope at early time, then turn around and roll down at late time. In 
other words, the exponentials exp(Pj • 0) decrease both in the past and in the future, at least if they 
explicitly appear in V, due to their associated integration constants Q being non-zero. Hence, for the 
moduli that are not completely frozen, the associated relative volumes of the blow-ups, e( TjJl ), always 
increases asymptotically at early and late times. We will confirm this general feature explicitly in 
our examples. 

5 A universal solution 

As a warm-up, we would first like to consider a universal solution where all bulk moduli and all 
blow-up moduli evolve in the same way. It appears that a natural way to obtain such a solution is 
to start with an Ansatz where the seven bulk moduli are equal, that is where there is only a single 
breathing mode. However, it turns out that such an Ansatz is incompatible with the equations of 
motion (4.10). The reason behind this is that not all bulk moduli are on the same footing because 
of their differing couplings to the blow-up moduli, as is evident from the coupling vectors in Table 1. 
Instead, we start with the slightly more general Ansatz 



where (j) is the breathing mode and ca are constants to be determined. For this Ansatz to be 
successful we need the seven equations of motion (4.10) for <pA to reduce to a single equation for <j>. 
This leads to the conditions 



4>a = ca4> 



(5.1) 



P(r,a) • c 



= const 



(5.2) 




(5.3) 



where we have defined 

$,a) = Yl C? T ,n,o) ' ( 5 - 4 ) 
n 

and is a constant. These conditions lead to a unique solution 

c = 1(4,4,4,3,3,3,3) (5.5) 

for the vector c in our Ansatz (5.1). It can further be shown that all conditions (5.2), (5.3) can be 
satisfied for specific choices of coefficients C(r,a) within a three-parameter family of solutions. The 
details of this are inessential for our subsequent discussion. For such a solution, the seven equations 
of motion (4.10) reduce to the single equation 

l0-2K 2 e- 2 ^ = O, (5.6) 

where fx is an arbitrary sign. This equation can be easily integrated once, and the integration constant 
can be fixed from the constraint (4.13). This leads to the first integral 

£ = |c| 2 Q</> 2 + C 2 e- 2 ^ (5.7) 

and the solution 

</» = ^lncosh(y) + |ln(M^) (5.8) 

y = ^i(r- Tl ) (5.9) 
l c l 

is a rescaled time variable and t\ is a constant. The seven bulk moduli are proportional to (/) and 
can be obtained by inserting this result into Eq. (5.1). From Eq. (4.4) one can then, in turn, obtain 
the solutions for the blow-up moduli. After another integration one finds 

Zi = -± (tanh(y) + 1) + z 0i (5.10) 



for the breathing mode <j>, where 



d t = ^ (5.11) 



where 

e 

and ZQi are independent integration constants. As advertised earlier, the blow-up moduli Z{ indeed 
approach constants asymptotically. More precisely, at early time, y — ► — oo, we have Z{ — > zoi and 
at late time, y — > oo, we have Zi — ► zu = z^i + di, that is, d{ is the distance by which the blow-up 
modulus Zi moves. We also note that the non-trivial evolution of Zi happens around the time y ~ 0. 
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Let us now interpret this solution. From Eq. (4.14) the orbifold volume is directly measured by 
the breathing mode 4> via 



We recall that our solution can only be trusted as long as q <C 1. Since, from Eq. (5.3), the pre-factor 
on the RHS of the above equation is of order one (for at least one i) we have to require the term in 
bracket be smaller than one. This leads to two cases, namely 

• €i <C 1 if |zoiAM 1 and In \zQi/di\ <C y <C — 1. In this time range the breathing mode evolves 
as (f) ~ — iiy + const. 

• €i -C 1 if <C 1 and 1 <C y <C In \di/z±i\. In this time range the breathing mode evolves 
as (f) ~ [iy + const. 

In both cases the breathing mode and, hence, the orbifold volume Vo, can increase or decrease 
depending on the choice of the sign However, independent of this choice, the relative blow-up 
volume £i will always leave the allowed range €j <C 1 when one of the limits of the given time ranges 
is approached. This means that after a finite proper time, both in the past and in the future, at least 
one of the blow-ups will take up a significant portion of the space and the approximation on which 
our Kahler potential (2.1) is based brakes down. In particular this means an evolution towards a 
state with small blow-ups at late time, y — > oo is not possible. 

6 Potential with a single exponential 

Let us now analyse the solutions more systematically, starting from simple patterns of evolution of 
the blow-up moduli and moving to more complicated ones. 

Certainly, the simplest possibility is to freeze all blow-up moduli by setting all Q = in Eq. (4.4). 
In this case, the effective potential (4.11) for the bulk moduli vanishes identically and the equations 
of motion (4.10) for <j)A can be easily integrated. In proper time t, related to r by Eq. (4.7), one 
easily finds as the general solution in this case 




(5.12) 



Further, from Eq. (4.15), the relative volume of a blow-up is given by 




(5.13) 




(6.1) 
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where to, T and Ua are constants. The expansion powers qA satisfy 

|q| 2 = | , (6.2) 

which follows from the constraint (4.13), and are otherwise arbitrary. These are simply solutions 
describing power-law evolution of the bulk moduli, which, in fact, are identical to the ones that can 
be obtained for M-theory on a seven-dimensional torus. 

We now move to the next more complicated case where blow-up moduli of only one particular 
type (r, a) evolve non-trivially and all the others have been set to constants. We write P = P( r ,o) 
and C 2 = 4 J2 n C(r,n,a) f° r simplicity of notation. The effective potential consists of only one term 
and takes the form 

V = C 2 e P ' . (6.3) 

This situation corresponds to an SU(2) Toda model so the general solution to Eq. (4.10) can be 
found. It is given by 

cj> = p« ln(ar) + (p (/) - p W ) ln(l + x s ) l ' s + k (6.4) 

subject to the constraints 

|p (i) | 2 = \ (6.5) 

5 = P-pW (6.6) 

P (/) = P (0 -^p^P (6-7) 

eX P( P " k ) = |p|2^2 " ( 6 - 8 ) 

x = ^ (6.9) 
is the rescaled proper time and to and T are constants. Note that Eqs. (6.5) and (6.7) imply that 

|P {/) | 2 = ^. (6-10) 

This solution is invariant under the exchange of pW and p^ and we remove this ambiguity by 
requiring S > in the positive-time branch and 8 < in the negative-time branch. The interpretation 
of these solutions is well-known [24, 25]. They interpolate between two, generally different, "free" 
solutions (6.1), one with q = pW at early time and one with q = p^ at late time. In these 
asymptotic regions the blow-up moduli Z( T , n ,a) are constant while they move around the time x ~ 1 
to facilitate the transition between the two free solutions. 
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Further, 



One can now compute the relative volume €( T ,a) °f the blow-ups by inserting the solution (6.4) 
into Eq. (4.15). As before, we find that the required condition t( T ,a) *C 1 is satisfied only in two finite 
time windows at x <C 1 and i> 1. At the endpoints of these windows the volume of the blow-ups 
becomes sizeable and control over the approximation is lost. 



7 More complicated cases 

An obvious generalisation of the previous example is to consider a situation that corresponds to an 
SU(2) n Toda model, for some integer n. This amounts to having a subset of blow-up moduli evolve 
non-trivially for which the associated characteristic vectors Pj are orthogonal, that is, 

P l • P, = 2% . (7.1) 

Inspection of Table 1 shows such cases can indeed be realised. For example, allowing the moduli 
of a certain type r to evolve non-trivially, while moduli of the other two types are being kept 
constant by setting Cr',n,a = in Eq. (4.4) for t' / r, leads to an SU(2) 3 Toda model with three 
exponentials appearing in the potential (4.11). From Table 1 there are a number of other options 
and in the following we will simply assume the existence of n vectors Pj satisfying (7.1) to cover all 
possibilities. 

The standard procedure to deal with such a system is to introduce a new constant basis = 
(ej, e a ) in field space, where i = 1, ■ ■ ■ , n and a = n + 1, ■ ■ ■ ,7 such that 

e, = -^fi (7.2) 

and the remaining vectors e a are chosen to complete the set to an orthonormal basis satisfying 

e A -e B = Sab ■ (7.3) 

The bulk fields can then be expanded as 

= PA*A , (7.4) 

A 

where pA is a new set of fields. The Lagrangian (4.12) and the Hamiltonian (4.13), written in terms 
of these new fields, take the form 

C = \Y.PA- y i H = \Y,Pa + V = E, (7.5) 

A A 
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with the potential 

y = 4^C, 2 e 2ft - (7.6) 

i 

Note that the modes p a have decoupled from this potential, which is, of course, one of the motivations 
behind introducing the basis ga- 

The general solution to the Lagrangian in (7.5) can be easily obtained as 

4> = J2pAeA (7.7) 



A 



with 



Pa = k a T + T a (7.8) 

Pi = v ^lncosh(y l )--^ln(^) , (7.9) 
where ta and Ua are constants and 

m = ^(r-n) (7.10) 

are rescaled time coordinates. The Hamiltonian constraint in (7.5) amounts to the condition 

Y,k\ = 4E. (7.11) 

A 

As expected the modes p a evolve freely. Each of the other modes pi evolves similarly to what has 
been found for the single SU{2) Toda model in the previous section. This can be seen explicitly 
by converting the above solution to proper time using (4.7) and comparing with Eq. (6.4). This 
means that each mode interpolates between two regions of simple power-law evolution at early and 
late time. The blow-up moduli Zi are constant in these asymptotic regions and their evolution at 
intermediate times facilitates the transition. This can be seen explicitly from their solution 

Zi = ^ (tanhfoi) + 1) + z 0i , (7.12) 

with zoi and di constants representing the value of the modulus at early time and its total change 
respectively. This solution is obtained by inserting the above solution for 4> into Eq. (4.4). 
From Eq. (4.15) one finds the relative volume 

of the blow-ups. As in all previous cases, can only be kept small, ej <C 1, for a finite proper time at 
either i/ 8 < -1 or > 1. For a valid solution, all e$ need to be small, which amounts to arranging 
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an overlap between those regions by choosing the time shifts appropriately. Outside this overlap 
region control of the approximation is lost as one or more blow-ups become large and take up a 
significant part of the internal space. 



There are also cases leading to a Toda model associated with a higher-rank simple group. Con- 
sider, for example, the three vectors 

P(a,i) = (-1,0,0,0,0,-1,0) (7.14) 
P m = (1,0,0,0,0,0,1) (7.15) 
P (7)3 ) = (0,-1,0,0,0,0,-1), (7.16) 

obtained from Table 1 with the sign choice M(a,i) = ~~ 1> M(/3, 2) = 1 an d ^(7,3) = ~~ 1 an d set all 
blow-up moduli with (r, a) different from the above to constants. The matrix A in Eq. (4.16) is then 
given by 

/ 2 -1 \ 

(A Ta )= -1 2-1 , (7.17) 

V - 1 2 / 

which is precisely the Cartan matrix of SU(4). Hence, we are dealing with an SU(4) Toda model. 

As before, the first step in solving this model is to introduce a new basis (e^) = (e r ,e„), where 
r = 1, 2, 3 and a = 4, . . . , 7, in field space. Here the three vectors e r are identified with (7.14)-(7.16) 
and the remaining four vectors e a are chosen to be orthonormal among themselves and orthogonal 
to all e r . Consequently, we have a basis satisfying 

e T ■ e a = A Tr7 , e T ■ e a = , e a ■ e b = 5 a b ■ ( 7 -18) 

Expanding the fields cf) 

(f, = Y J PAe A (7.19) 

A 

as before one finds for the Lagrangian (4.12) 



c = \ E a t,°ptp* + -JZpI- v ( 7 - 2 °) 

T,cr a 

with potential 

y = 4 E^ 2ex pfe j4 -^ • ( 7 - 21 ) 

The modes p a are decoupled from this potential and their equations of motion immediately lead to 
the general solution 

Pa = k a (r-T a ) , (7.22) 
14 



where k a and r a are constants. Following the methods of Ref. [21], the solutions for the other modes 
are obtained as 

e'Pr = b T (X) exp (A • (kr - r)) , (7.23) 
AeA T 

where k and r are constant vectors. The three sets A T contain the weights of the fundamental 
representation 4, the vector representation 6 and the anti- fundamental representation 4 of SU(4), 
respectively, and are explicitly given by 

Ai = {(100), (-110), (0 - 11), (00 - 1)} (7.24) 
A 2 = {(010), (1-11), (-101), (10-1), (-11-1), (0-10)} (7.25) 
A 3 = {(001), (01 - 1), (1 - 10), (-100)} . (7.26) 

Finally, b T (X) represents a set of constants which depends on £ T as well as on k. They can be 
calculated by inserting the solutions (7.23) into the equations of motion but their explicit form will 
not be of any relevance here. 

The important feature of (7.23) is that asymptotically one of the exponentials in the sum will 
dominate leading to a power-law evolution. In these regions, the blow-up moduli are approximately 
constant and the relative blow-up volumes increase until they approach values of order one where the 
approximation underlying our analysis breaks down. Hence, this more complicated Toda solution 
also conforms with our general expectation of asymptotically increasing blow-up volumes. 
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